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ABSTRACT 

The coming few years are likely to witness a dramatic increase in high quality Sn 
data as current surveys add more high redshift supernovae to their inventory and 
as newer and deeper supernova experiments become operational. Given the current 
variety in dark energy models and the expected improvement in observational data, 
an accurate and versatile diagnostic of dark energy is the need of the hour. This 
paper examines the Statefinder diagnostic in the light of the proposed SNAP satellite 
which is expected to observe about 2000 supernovae per year. We show that the 
Statefinder is versatile enough to differentiate between dark energy models as varied 
as the cosmological constant on the one hand, and quintessence, the Chaplygin gas and 
braneworld models, on the other. Using SNAP data, the Statefinder can distinguish 
a cosmological constant (w = —1) from quintessence models with w —0.9 and 
Chaplygin gas models with k ^ 15 at the 3<r level if the value of Q m is known exactly. 
The Statefinder gives reasonable results even when the value of fi m is known to only 
~ 20% accuracy. In this case, marginalizing over f2 m and assuming a fiducial LCDM 
model allows us to rule out quintessence with w ^ —0.85 and the Chaplygin gas 
with k ^ 7 (both at 3cr). These constraints can be made even tighter if we use the 
Statefmders in conjunction with the deceleration parameter. The Statefinder is very 
sensitive to the total pressure exerted by all forms of matter and radiation in the 
universe. It can therefore differentiate between dark energy models at moderately 
high redshifts of z < 10. 



Key words: cosmology: theory — cosmological parameters — statistics 



1 INTRODUCTION 



Supernova observations (Riess et al. 1998; Perlmutter et al. 
1999), when combined with those of the cosmic microwave 
background (Benoit et al. 2003) and gravitational cluster- 
ing (Percival 2002), suggest that our Universe is (approxi- 
mately) spatially flat and that an exotic form of negative- 
pressure matter called 'dark energy' (DE) causes it to ac- 
celerate by contributing as much as 2/3 to the closure den- 
sity of the universe - the remaining third consisting of non- 
relativistic dark matter and baryons. The simplest example 
of dark energy is the cosmological constant (A), with asso- 
ciated mass density 



p A = 6.44xl0- 30 (^)(A) 2 gcm -3 ; (1) 

where h is the Hubble constant Ho in terms of 100 km s _1 
Mpc -1 andfi A = 0.7±0.1, /i = 0.7±0.1. Although the cold 
dark matter model with a cosmological constant (hereafter 
LCDM) provides an excellent explanation for the accelera- 
tion phenomenon and other existing observational data, it 
remains entirely plausible that the dark energy density is 
weakly time dependent (see the reviews Sahni & Starobin- 
sky 2000; Peebles & Ratra 2003). Moreover, it is natural to 
suggest (in complete analogy with what has been done in the 
case of another type of 'dark energy' responsible for driving 
the expansion of the universe during an inflationary stage in 
the early universe) that the dark energy which we observe 
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today might really be dynamical in nature and origin. This 
means that a completely new form of matter is responsible 
for giving rise to the second inflationary regime which we 
are entering now. 

Many models of dark energy have been proposed; in 
fact, any inflationary model (even a 'bad' one, i.e. with- 
out a 'graceful exit' to the subsequent radiation-dominated 
Friedmann- Robertson- Walker (FRW) stage) may be used for 
this purpose if one assumes different values for its micro- 
scopic parameters. The simplest of these models rely on a 
scalar field minimally interacting with Einstein gravity - 
quintessence (Ratra & Peebles 1988; Peebles & Ratra 1988; 
Frieman et al. 1995; Caldwell, Dave & Steinhardt 1998), 
and bear an obvious similarity with the simplest variants 
of the inflationary scenario. Inclusion of a non-minimal cou- 
pling to gravity in these models together with further gen- 
eralization leads to models of dark energy in a scalar-tensor 
theory of gravity (see Boisseau et al. 2000, and references 
therein). Other models invoke matter with unusual proper- 
ties such as the Chaplygin gas (Kamenshchik, Moschella & 
Pasquier 2001) or k-essence (Armendariz-Picon, Mukhanov 
& Steinhardt 2000). Still others generate cosmic accelera- 
tion through topological defects (Bucher & Spergel 1999) 
or quantum vacuum polarization and particle production 
(Sahni & Habib 1998; Parker & Raval 1999). Lately it has 
been noticed that higher dimensional 'braneworld' models 
could account for a late-time accelerating phase even in 
the absence of matter violating the strong energy condi- 
tion (Dvali, Gabadadze & Porrati 2000; Deffayet, Dvali & 
Gabadadze 2002; Deffayet et al. 2002; Sahni & Shtanov 
2002; Alam & Sahni 2002) (see Sahni 2002, for a recent re- 
view of dark energy models) . It is especially interesting that 
in the latter class of models 'dark energy' need not be an 
energy of some form of matter at all, but can have an en- 
tirely geometrical origin. Moreover, in these models the basic 
gravitational field equations do not have the Einstein form 

Ra/3 — -g a fjR = 8nG (T al 3\ +T a fj\ +T a fj\ )(2) 

2 v 1 matter ' \ radiation 1 \DE> 

(c = h — 1 is assumed here and below), and therefore 
the notions of 'energy density' and 'pressure' of DE loose 
their exact fundamental sense and become ambiguous and 
convention-dependent. A major ambiguity arises in mod- 
els of scalar-tensor gravity as well as in braneworld models 
both of which contain interaction terms between dark energy 
and non-relativistic matter. Interpreting such models within 
the Einstein framework (2) leads to the following dilemma: 
should these interaction terms be ascribed to dark matter 
(hence to T a g\ ma tter in (2)) or to dark energy (to T aj )\ DE ) 
? Our answer to this question has the potential to alter the 
properties of dark energy including its density and pressure 
and hence also its equation of state. In marked contrast to 
such ambiguities which could arise if we are not careful with 
our usage of the term 'equation of state', the expansion fac- 
tor of the universe in the physical frame a(t), when expressed 
through the Hubble parameter H = -, is an unambiguous, 
fundamental and readily measurable quantity. 

Given the rapidly improving quality of observational 
data and also the abundance of different theoretical mod- 
els of dark energy, the need of the hour clearly is a robust 
and sensitive statistic which can succeed in differentiating 
cosmological models with various kinds of dark energy both 



from each other and, even more importantly, from an exact 
cosmological constant. In view of the non-fundamentality of 
the notions of DE density and pressure pointed out above, 
we prefer to work with purely geometric quantities. Then 
such a sensitive diagnostic of the present acceleration epoch 
and of dark energy could be the statefinder pair {r, s}, re- 
cently introduced in Sahni et al. (2003). The statefinder 
probes the expansion dynamics of the universe through 
higher derivatives of the expansion factor a & a . Its im- 
portant property is that {r,s} — {1,0} is a fixed point for 
the flat LCDM FRW cosmological model. Departure of a 
given DE model from this fixed point is a good way of es- 
tablishing the 'distance' of this model from flat LCDM. As 
we will show in this paper, the statefinder successfully differ- 
entiates between rival DE models and, when combined with 
SNAP supernova data, can serve as a versatile and powerful 
diagnostic of dark energy. 

The paper is organized as follows. In the next section 
we briefly review some theoretical models of dark energy. 
The behaviour of the statefinder pair for these models is 
discussed in Section III while the nature of data expected to 
become available from the SNAP experiment is the subject 
of Section IV. Section IV also discusses model-independent 
parametric reconstructions of dark energy. Our conclusions 
are presented in section V. 

2 DARK ENERGY MODELS AND THE 
ACCELERATION OF THE UNIVERSE 

The rate of expansion of a FRW universe and its acceleration 
are described by the pair of equations 

„„ StyG k 

H = —2^1*-*' 

i 

a 4:7iG , „ x 

- = 3~2>i + 3pi)> (3) 

i 

where the summation is over all matter fields contributing 
to the dynamics of the universe. Clearly a necessary (but 
not sufficient) condition for acceleration (o > 0) is that at 
least one of the matter fields in (3) violate the strong en- 
ergy condition p + 3p > 0. If for simplicity we assume that 
the dark energy pressure and density are related by the sim- 
ple linear relation p = wp, then w < —1/3 is a necessary 
condition for the universe to accelerate. The acceleration of 
the universe can be quantified through a dimensionless cos- 
mological function known as the 'deceleration parameter' 
q = —a/aH 2 , equivalently 

q(x) = ^x-l, x = l + z, (4) 

where q < describes an accelerating universe, whereas q ^ 
for a universe which is either decelerating or expanding 
at the 'coasting' rate a oc t. As it will soon be shown, the 
deceleration parameter on its own does not characterize the 
current accelerating phase uniquely. The presence of a fairly 
large degeneracy in q(z) is reflected in the fact that rival 
dark energy models can give rise to one and the same value 
of go at the present time. This degeneracy is easily broken 
if, as demonstrated in section 3, one combines q(z) with 
one of the statefinders r(z), s(z). The diagnostic pairs {r, q} 
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and {s, q} provide a very comprehensive description of the 
dynamics of the universe and consequently of the nature of 
dark energy. 

Now let us come to the issue of defining the energy 
density and pressure of DE. In view of the ambiguities dis- 
cussed in the Introduction, we shall define pde and pdb by 
making use of the Einstein interpretation of gravitational 
field equations (not to be confused with the notion of the 
Einstein frame which is used in scalar-tensor and string the- 
ories of gravity!). Namely, we assume that the gravitational 
field equations in a single- metric theory of 3+1 gravity can 
be formally written in the form (2) where the Einstein ten- 
sor standing in the left-hand side is defined with respect to 
the physical space-time metric. All other terms are trans- 
ferred to the right-hand side. Next, we subtract the energy- 
momentum tensor of dust (CDM + baryons) from the total 
energy-momentum tensor of matter and call the remaining 
part 'the effective energy-momentum tensor of dark energy' 
(in the Einstein interpretation) . Combining this prescription 
with Eq. (3) and in the absence of spatial curvature, the en- 
ergy density and pressure of dark energy can be denned as: 

3H 2 

Pde = 



Pde 



pcritical — Pm — g^Q (1 ~ I 

H 2 , F 



(5) 



where p cr iticai = 3H 2 /8ttG is the critical density associated 
with a FRW universe. An important consequence of using 
this approach is that the ratio w DE = Pde/pde can be 
expressed in terms of the deceleration parameter 



Wed (x) 



2q(x) 



3(i-n m (x)) 



(2x/3) dlnH / dx-1 
1 - (H /H) 2 n m x 3 



(6) 



Following the above prescription we get standard results 
for the cosmological constant and quintessence (for instance 
we recover Eq. (9)). However the same cannot be said of 
braneworld models since the Hubble parameter for the lat- 
ter contains interaction terms between matter and dark en- 
ergy (see for instance Eqs. (22), (23)) and therefore does not 
subscribe to the Einsteinian format (2) & (3). One can how- 
ever extend the above definition of w e g to non-Einsteinian 
theories by defining dark energy density to be the remain- 
der term after one subtracts the matter density from the 
critical density in the Einstein equations. It should be em- 
phasised that, according to this prescription all interaction 
terms between matter and dark energy (such terms arise in 
scalar-tensor and brane models) are attributed solely to dark 
energy. Therefore u> e ff (z) defined according to (6) is an effec- 
tive equation of state in these models and not a fundamental 
physical entity (as it is in LCDM, for instance). 

In this connection we should also stress that the prop- 
agation velocity of small inhomo geneities in d ark energy is 
generically neither ^wde, nor U dpoE/dpoE- Therefore al- 
though w(z) is an important physical quantity it does not 
provide us with an exhaustive description of dark energy and 
its use as a diagnostic should be treated with some caution. 
(In this paper we will restrict ourselves to a spatially flat 
FRW model and will not consider inhomogeneous perturba- 
tions on this background.) 

We now highlight a few popular candidates for dark 
energy which shall be the focus of our discussion in this 
paper. 



• Cosmological Constant. Perhaps the simplest model 
for dark energy is a cosmological constant A, whose energy 
density remains constant with time pa = A/8ttG = — pa, 
and which has an equation of state wa = — 1. A universe 
consisting of matter in the form of dust and the cosmo- 
logical constant is popularly known as LCDM, the Hubble 
parameter for this model has the form 



H(z) = H [O m (l + zf + 1 - 0, 



I 1/2 



(7) 



• Quiessence. The next simplest form of dark energy 
after the cosmological constant is provided by models for 
which the equation of state is a constant w = constant ^ — 1. 
For this form of dark energy, which we call 'quiessence' 



w) 



1/2 



(8) 



H{z) = Ho [O m (l + zf + Q X {1 + z) 3(1+ 

For w — — 1 we recover the limiting form (7). Important 
examples of quiessence include a network of non-interacting 
cosmic strings (w = —1/3) and domain walls (w = —2/3). 
Quiessence in a FRW universe can also be produced by a 
scalar field (quintessence) which has the potential V{<j>) oc 

-2(l + m) 

sinh «i (C(j> + D), with appropriately chosen values of 
C and D (see Sahni & Starobinsky 2000; Urena-Lopez & 
Matos 2000). 

Usually the dark energy equation of state depends upon 
time. We call such more generic models kinessence. 

• Quintessence The simplest example of kinessence is 
provided by quintessence - a self-interacting scalar field 
which couples minimally to gravity. Its density, pressure and 
equation of state are given by 



P4> 



U 2 + v(4>) 



2 

> -1 



(9) 



Scalar field evolution is governed by the equation of motion 

+ 3#0+^=O, (10) 
dg> 

where 

s.^r: r i .„ i 

(11) 



R 2 _ 8ttG 



Pom(i + zy + -<f> 2 + v(<i>) 



It is clear from (9) that w < -1/3 provided <j> 2 < V(<f>). Mod- 
els with this property can lead to an accelerating universe 
at late times. An important subclass of quintessence models 
displays the so-called 'tracker' behaviour during which the 
ratio of the scalar field energy density to that of the mat- 
ter/radiation background changes very slowly over a sub- 
stantial period of time. Models belonging to this class sat- 
isfy V"V/(V') 2 > 1 and approach a common evolutionary 
'tracker path' from a wide range of initial conditions. As 
a result, the present value of dark energy in tracker mod- 
els is to a large extent (though not entirely) independent of 
initial conditions and is determined by parameters residing 
only in its potential - as in the case of the cosmological con- 
stant (for a brief review of tracker models see Sahni 2002). 
In this paper we will focus our attention on the tracker po- 
tential V{4>) oc <fT a , a ^ 1 which was originally proposed 
in Ratra & Peebles (1988). For this potential, the region 
of initial conditions for cf> f° r which the tracker regime has 
been reached before the end of the matter-dominated stage 
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is 4>in <S Mp = 1/ vG, and the present value of quintessence 
is 4>(t ) ~ M P . 

For all quintessence models w ^ — 1, and this inequality is 
saturated only if <j> = dV/d<j> = 0. In order to obtain w < — 1 
matter must violate the strong energy condition p + 3p ^ 0, 
for some duration of time. It should be noted that DE with 
w < — 1 is not excluded by observations (see Melchiorri et al. 
2002, for a recent investigation). However in order to have 
w < — 1 one must look beyond quintessence models. Models 
based on scalar-tensor gravity (Boisseau et al. 2000) can 
have w < — 1, so too can braneworld models (see Sahni & 
Shtanov 2002 for a discussion of this issue and Alam & Sahni 
2002 for a comparison of braneworld models with observa- 
tional data). 

• Chaplygin gas. An interesting alternate form of dark 
energy is provided by the Chaplygin gas (Kamenshchik et al. 
2001; Bilic, Tupper & Viollier 2002; Fabris, Goncalves & 
de Souza 2002; Gorini, Kamenshchik & Moschella 2003; Al- 
caniz, Jain & Dev 2003; Avelino et al. 2003) which obeys 
the equation of state 



Pc 



-A/p c 



(12) 



The energy density of the Chaplygin gas evolves according 
to 



Pc = ^JA + B(l + z) e 



(13) 



from where we see that p c — > \/A as z — > —1 and p c — > 
\/B(l + z) 3 as z > 1. Thus, the Chaplygin gas behaves 
like pressureless dust at early times and like a cosmological 
constant during very late times. Note, however, that Chap- 
lygin gas at z >• 1 is not simply a new kind of CD M if we 
examine its inhomogeneities ( i. e. if we apply this hydrody- 
namical equation of state to the inhomogeneous case, too)! 
In contrast to CDM a nd baryo ns, the sound velocity in the 
Chaplygin gas v c = U dp c /dp c — \fA/ p c quickly grows oc t 2 
during the matter-dominated stage and becomes of the or- 
der of the velocity of light at present (it approaches light 
velocity asymptotically in the distant future ). Thus, from 
the point of view of inhomogeneities, the properties of the 
Chaplygin gas during the matter-dominated epoch are very 
unusual and resemble those of hot dark matter which has a 
large Jeans length, despite the fact that the Chaplygin gas 
formally carries negative pressure. 

The Hubble parameter for a universe containing cold dark 
matter and the Chaplygin gas is given by 



H(z) 



Ho 



n m (i+^) 3 +— \/4+(i+^) 6 

Hi y JD 



1/2 



(14) 



where k = pom/V~B. It is easy to see from (14) that 
K= ^1( Z _ >00 ). (15) 

Pc 

Thus, k defines the ratio between CDM and the Chaplygin 
gas energy densities at the commencement of the matter- 
dominated stage. It is easy to show that 



(16) 



In the limiting case when A = 0, the Chaplygin gas becomes 
indistinguishable from dust-like matter (if we examine its be- 



haviour in an unperturbed FRW background) . This limiting 
case corresponds to 



^T7 



(17) 



and is shown as the outer envelope (dashed) to the Chap- 
lygin gas models in Figures la,b. In the other limiting case 
B = 0, the Chaplygin gas reduces to the cosmological con- 
stant. 

The fact that the sound velocity in the Chaplygin gas 
is not small during the matter-dominated stage and be- 
comes very large towards its end suggests that the parame- 
ter k should be large in order to avoid damping of adiabatic 
perturbations. This requires A ^> B. Recent investigations 
which look at Chaplygin gas models in the light of galaxy 
clustering data and CMB anisotropies show that this ob- 
servation is correct if the equation of state p c oc — l/p c is 
assumed to be universally valid (Carturan & Finelli 2002; 
Sandvik et al. 2002; Bean & Dore 2003). In our paper we 
consider the Chaplygin gas equation of state to be a phe- 
nomenological description of dark energy in a FRW back- 
ground and do not assume that it remains true for pertur- 
bations. However, the fact that k should be large for vi- 
able models will appear in our results, too. Finally let us 
point out that the Chaplygin gas may be considered to be 
a specific case of k-essence with a constant potential and 
the Born-Infeld kinetic term. To illustrate this consider the 
Born-Infeld lagrangian density 

£ = -Vo V / l-4>,^ , (18) 

where <f>^ = d(j>/dx' 1 . For time-like one can define a four 
velocity 



(19) 



this leads to the standard form for the hydrodynamical 
energy-momentum tensor 



Ty. v = (p + p)u^u„ — pg^ v , 

where (Frolov, Kofman & Starobinsky 2002) 



P = 



V 



, p = -Vos/l 



(20) 



(21) 



and we find that we have recovered (12) with A = V 2 . 

• Braneworld models. Braneworld models provide an 
interesting alternative to dark energy model building. Ac- 
cording to this higher dimensional world view, we live on 
a 3+1 dimensional brane ('brane' being a multidimensional 
generalization of 'membrane') which is either embedded in 
or bounds a higher dimensional space-time. The simplest 
example of a braneworld which can lead to late-time accel- 
eration is the model suggested by Deffayet et al. (2002) (we 
shall henceforth refer to this model as the DDG model). 



H 



8tyGp„ 



1 



(22) 



where l c — m 2 /M 3 is a new length scale and m and M 
refer respectively to the four and five dimensional Planck 
mass (l c = 2r c in the terminology of Deffayet et al. 2002). 
The acceleration of the universe in this model is not caused 
by the presence of 'dark energy' but due to the fact that 
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general relativity is formulated in 5 dimensions instead of 
the usual 4. One consequence of this is that gravity becomes 
five dimensional on length scales R > l c = 2H^ 1 (1 — f2 m ) _1 . 
A more general class of braneworld models is described by 
(Sahni & Shtanov 2002) 

= O m (l+z) 3 + CT + 20 i r 



y/n^jl + z) 3 + » CT + + Q Ab , (23) 

where Ab is the bulk cosmological constant, a is the brane 
tension and 



POm 

3m 2 11$ 



3m 2 # 2 ' 



A b 



(24) 



It is easy to see that l c can be of the same order as the 
Hubble radius l c ~ H' 1 if M ~ 100 MeV. On short length 
scales r <C l c and at early times, one recovers general rel- 
ativity, whereas on large length scales r ^> l c and at late 
times brane-related effects begin to play an important role. 
It is interesting that brane-inspired effects can lead to the 
late time acceleration of the universe even in the complete 
absence of a matter source which violates the strong energy 
condition p + 3p ^ (Deffayet et al. 2002; Sahni & Shtanov 
2002). 

The dimensionless value of the brane tension Q<? is deter- 
mined by the constraint relation 



n m + a, =f 2VWi y / i-n K + n Ah = 1. 



(25) 



The underlined terms in (23) & (25) make braneworld 
models different from standard FRW cosmology. Indeed by 
setting fij = (23) reduces to the LCDM model 



On(l + z) 3 



(26) 



which describes a universe containing matter and a cosmo- 
logical constant (7). The two signs in (23) correspond to the 
two separate ways in which the brane can be embedded in 
the higher dimensional bulk. As shown in Sahni & Shtanov 
(2002), taking the upper sign in (23) and (25) leads to the 
model called BRANE1, while the lower sign in (23) and (25) 
results in BRANE2. 

Three important classes of braneworld models deserve 
special mention: 

(i) BRANE1 models have an effective equation of state 
which is more negative than that of the cosmological con- 
stant w < — 1. 

(ii) BRANE2 models have w ^ — 1. For parameter values 
Q<? — QA b = 0, BRANE2 coincides with the dark energy 
model discussed in Eq. (22). 

(iii) A class of braneworld models, called 'disappearing 
dark energy' (DDE) (Sahni & Shtanov 2002; Alam & 
Sahni 2002), have the important property that the current 
acceleration of the universe is a transient phase which is 
sandwiched between two matter dominated epochs. These 
models do not have horizons and therefore help to rec- 
oncile an accelerating universe with the demands of the 
string/M-theory (Sahni 2002) (as well as any theory which 
requires dark energy to decay in the future and transform 
into matter with w > —1/3). 



Finally we note that, for a spatially flat universe, the lu- 
minosity distance for all models discussed above is given by 
the simple expression 



Dl(z) 
1 + z 



dz' 
H(z') 



(27) 



where H(z) is given by (7) for LCDM, by (8) for quiessen- 
cence, by (11) for quintessence, by (14) for the Chaplygin 
gas and by (23) for the braneworld models. 



3 THE STATEFINDER DIAGNOSTIC 

As we have seen above, dark energy has properties which can 
be very model dependent. In order to be able to differenti- 
ate between the very distinct and competing cosmological 
scenarios involving dark energy, a sensitive and robust di- 
agnostic (of dark energy) is a must. Although the rate of 
acceleration/deceleration of the universe can be described 
by the single parameter q — —ii/aH 2 , a more sensitive dis- 
criminator of the expansion rate and hence dark energy can 
be constructed by considering the general form for the ex- 
pansion factor of the Universe 



a(t) = a(t ) 



+a\ Q (t 



-to)+-£(t- 



to) 



a 



tof + . 



.(28) 



In general, dark energy models such as quiessence, 
quintessence, k-essence, braneworld models, Chaplygin gas 
etc. give rise to families of curves a(t) having vastly differ- 
ent properties. Since we know that the acceleration of the 
universe is a fairly recent phenomenon (Benitez 2002; Riess 
2001; Sahni & Starobinsky 2000) we can, in principle, con- 
fine our attention to small values of \t — to\ in (28). We have 
shown in Sahni et al. (2003) that a new diagnostic of dark 
energy called statefinder can be constructed using both the 
second and third derivatives of the expansion factor. The 
second derivative is encoded in the deceleration parameter 
which has the following form in a spatially fiat universe: 



a 

"alP 



(i + 3wn x ) 



(29) 



The statefinder pair {r, s}, defines two new cosmological pa- 
rameters (in addition to H and q): 



a 9w_ „ . 3 w 



r-l 



3(9-1/2) 



1 w 

1 + w 's1Fh 



(30) 
(31) 



Clearly an important requirement of any diagnostic is 
that it permits us to differentiate between a given dark en- 
ergy model and the simplest of all models - the cosmolog- 
ical constant A. The statefinder does exactly this. For the 
LCDM model, the value of the first statefinder stays pegged 
at r = 1 even as the matter density evolves from a large 
initial value (Q m — 1, t <JC to) to a small late-time value 
(f2 m — » 0, t ^> to). It is easy to show that {r, s} = {1, 0} is 
a fixed point for LCDM. 

The second statefinder s has properties which comple- 
ment those of the first. Since s does not explicitly depend 
upon either Qx or fi m , many of the degeneracies which 
are present in r are broken in the combined statefinder 
pair {r, s}. For models with a constant equation of state 
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Figure 1. The left panel (a) shows the time evolution of the statefinder pair {r,s} for quintessence models and the Chaplygin gas. 
Quintessence models lie to the right of the LCDM fixed point (r = l,s = 0) (solid lines represent tracker potentials V = Vo/<j> a , 
dot-dashed lines representing quicsscncc with constant equation of state w). For quiessence models, s remains constant at 1 + w while 
r declines asymptotically to 1 + §«>(1 + w). For tracker models, s monotonically decreases to zero whereas r first decreases from unity 
to a minimum value, then rises to unity. These models tend to approach the LCDM fixed point (r = 1, s = 0) from the right at t — * oo. 
Chaplygin gas models (solid lines) lie to the left of the LCDM fixed point. For these models, k is the ratio between matter density and 
the density of the Chaplygin gas at early times. For all Chaplygin gas models, s monotonically increases to zero from -1, whereas r first 
increases from unity to a maximum value, then decreases (to unity) . The dashed curve in the lower right is the envelope of all quintessence 
models, while the dashed curve in the upper left is the envelope of Chaplygin gas models (the latter is described by k = f2 m /l — f2 m ). 
The region outside the dashed curves is forbidden for both classes of dark energy models. The right panel (b) shows the time evolution 
of the pair {r, q}, where q is the deceleration parameter. It is important to note that the the solid line, which corresponds to the time 
evolution of the LCDM model, divides the r — q plane into two halves. The upper half is occupied by Chaplygin gas models, while the 
lower half contains quintessence models. All models diverge at the same point in the past (r = 1, q = 0.5) which corresponds to a matter 
dominated universe (SCDM), and converge to the same point in the future (r = 1, q = — 1) which corresponds to the steady state model 
(SS) — the de Sitter expansion (LCDM — » SS as t — » oo and Q m — * 0). The dark dots on the curves show current values {rrj,so} (left) 
and {ro,<?o} (right) for different dark energy models. In all models, f2 m = 0.3 at the current epoch. In both panels quicsscncc is shown 
as dot-dashed while dashed lines mark envelopes for Chaplygin gas (upper) and quintessence (lower). 



(quiessence) s — 1 + w = constant, while the statefinder r is 
time-varying. For models with time-dependent equation of 
state (kinessence) , both r and s vary with time. As we will 
show in this paper, the statefinder pair {r, s} can easily dis- 
tinguish between LCDM, quiessence and kinessence models. 
It can also distinguish between more elaborate models of 
dark energy such as braneworld models and the Chaplygin 
gas (see also Sahni et al. 2003, Gorini et al. 2002). Inter- 
estingly, as demonstrated in section 5, the statefinder pair 
{s, q} proves to be an even better diagnostic of dark energy 
than {r, s}. 

The statefinders r and s can be easily expressed in terms 
of the Hubble parameter H(z) and its derivatives as follows: 




s{x) = 



3(g(x) - 1/2) ' 



(32) 



where x = 1+z and H is given by (7), (8), (11), (14), (23) for 
the different dark energy models discussed in the previous 
section. 

In figure 1(a) we show the time evolution of the 
statefinder pair {r, s}. We find that the vertical line at s = 
effectively divides the r — s plane into two halves. The left 
half contains Chaplygin gas (CG) models which commence 
their evolution from r = 1, s = —1 and end it at the LCDM 
fixed point (r = l,s = 0) in the future. The quintessence 
models occupy the right half of the r — s plane. These models 
commence their evolution from the right of the LCDM fixed 
point and, like CG, are also attracted towards the LCDM 
fixed point in the future. For quiessence models, r decreases 
monotonically to 1 + §to(l + w) while s remains constant 
at s = 1 + w. For kinessence models, on the other hand, s 
decreases monotonically to zero, while r first decreases to a 
minimum value then increases to unity. The region below the 
curve r = 1 + |s(s — 1) is disallowed for quintessence models 
whereas CG models with k < fi m /(l — 0, m ) are excluded. 
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Figure 2. Trajectories in the statefinder plane {r,q} for the 
brancworld models discussed in (23). BRANE1 models have 
w ^ — 1 gcnerically, whereas BRANE2 models have w ^ — 1. 
The closed loop represents DDE in which the acceleration of 
the universe is a transient phenomenon. For braneworld mod- 
els, parameter values are as follows. BRANE1: solid curves above 
LCDM; top to bottom: Q m = 0.6, Q t = 6.0, n m = 0.5, S7; = 2.0, 
n m = 0A,Q l = 0.5. BRANE2: solid curves below LCDM; 
top to bottom: tt m = 0.3,^ = 0.05, n m = 0.2, = 0.25, 
tt m = 0.1, n t = 0.45. The thick solid curve in BRANE2 corre- 
sponds to the DDG model discussed in (22) with f! m = 0.24. 
For BRANE1 and BRANE2 models U Ab = (i.e. there is no 
cosmological constant in the bulk.) Dark dots indicate the cur- 
rent value of {r, q} for the models. All models are in reasonable 
agreement with current supernova data. For DDE, from outer 
to inner loops, f2 m = 0.05, QA b = 4.9, Q m = 0.15, U\ b = 1.4, 

n m = 0.20, n Ab = i.i. 



It is interesting that the second statefinder, s, is positive for 
quintessence models, but negative for the CG. Similarly the 
first statefinder, r, is < 1 (> 1) for quintessence (CG). The 
distinctive trajectories which quiessence, quintessence and 
CG follow in the r — s plane demonstrates quite strikingly 
the contrasting behaviour of dark energy models. 

The separation between distinct families of dark en- 
ergy models is also very pronounced when we analyze evo- 
lutionary trajectories using the statefinder pair {r, g} shown 
in Fig. 1(b) and Fig. 2. Fig. 1(b) shows the evolution of 
quintessence and CG models in r — q space, while Fig. 2 
shows the evolution of the braneworld models discussed in 
(23). In Fig. 1(b) the LCDM model effectively divides the 
r — q space in half, separating quintessence models (bottom- 
half) from the Chaplygin gas (top-half). From this figure we 
clearly see that all dark energy models commence evolv- 
ing from the same point in the past (g = 0.5, r — 1), 
which corresponds to a matter dominated SCDM universe. 
Quintessence, LCDM and the Chaplygin gas all end their 
evolution at the same common point in the future (q = 



— l,r — 1), which corresponds to steady state cosmology 
(SS) - the de Sitter expansion. In Fig. 2 the LCDM model 
separates BRANE1 models (which have w e g < — 1) from 
BRANE2 models as well as DDE models. BRANE2 models 
have w e ff ^ — 1 generically, whereas DDE models consist of 
a transient accelerating regime which is sandwiched between 
two matter dominated epochs. Thus DDE both begins and 
ends its evolution at the SCDM point {r, q} = {1,0.5} and 
its r — q space trajectory is a loop ! BRANE1 and BRANE2 
models on the other hand, commence evolving at the SCDM 
point and tend to SS in the future. Fig. 1(b) and Fig. 2 
clearly demonstrate that the deceleration parameter can- 
not on its own differentiate between rival models of dark 
energy. The degeneracy which afflicts q(z) clearly also af- 
flicts the equation of state w(z), since both q and w are 
related through (6). We therefore feel we have convincingly 
demonstrated that the statefinders can successfully differen- 
tiate between competing dark energy models as diverse as 
LCDM, quintessence, braneworld models and the Chaplygin 
gas. Statefinders can also be applied to other interesting can- 
didates for dark energy including bigravity models (Damour, 
Kogan & Papazoglou 2002), generalized Chaplygin gas (Ka- 
menshchik, Moschella & Pasquier 2001; Bento, Bertolami 
& Sen 2002), k-essence (Armendariz-Picon, Mukhanov & 
Steinhardt 2000) scalar-tensor theories etc. 

Finally we draw the readers attention to the following 
elegant relationship which exists between the statefinders 
on the one hand, and the total density p = £\ pi and total 
pressure p — ^\ pi in the universe: 

1 |p r _ 1= Hp+p)i s= (p±p)£ (33) 

2 2p ' 2p p ' p p y ' 

From Eq. (33) we see that the statefinder s is exceedingly 
sensitive to the total pressure p. This has some interesting 
consequences. At early times the presence of radiation en- 
sures that the total pressure in the universe is positive. Much 
later, the universe begins to accelerate driven by the nega- 
tive pressure of dark energy. In between these two asymp- 
totic regimes, deep in the matter dominated epoch, a stage 
is reached when the (negative) pressure of dark energy is 
exactly balanced by the positive pressure of radiation. At 
this precise moment of time p ~ and s — * oo ! For LCDM 
this pressure balance is achieved at 2, ~ 10, consequently 
\s\ 1 when 2 ~ 2*. It can be shown that the redshift 2* 
(at which p = 0) is quite sensitive to the form of dark en- 
ergy. We therefore find that the statefinder s diagnoses the 
presence of dark energy even at high redshifts when the con- 
tribution of DE to the total energy budget of the universe 
is insignificant ! 



4 MODEL INDEPENDENT 

RECONSTRUCTION OF COSMOLOGICAL 
PARAMETERS FROM SNAP DATA 

4.1 The cosmological reconstruction of dark 
energy properties 

Cosmological reconstruction is an effective statistical tech- 
nique which can be used in situations where a large number 
of theoretical models are to be compared with observations. 
Instead of estimating relevant parameters for each model 
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separately, we can choose a model-independent fitting func- 
tion and perform a maximum likelihood parameter estima- 
tion for it. The resultant confidence levels can be used to rule 
out or accept the different models available. This technique 
is effective here because, as discussed in Section 2, a wide 
range of theoretical models have been suggested to explain 
dark energy. 

The basis of cosmological reconstruction rests in the 
observation that the expression for the luminosity distance 
(27) can be easily inverted (Starobinsky 1998; Huterer & 
Turner 1999; Nakamura & Chiba 1999): 



H(z) 



d_ 

dz 



Dl(z) 
l + z 



(34) 



Thus, from mathematical point of view, any given Dl(z) de- 
fines H(z). Eqs. (5) and (6) can then be used to obtain the 
dark energy density and the associated equation of state. 
Similarly the statefinder pair {r, s} can be determined by 
employing Eq (32) together with Eq (29) . However, in prac- 
tice the derivative with respect to z may not be simply 
performed since Dl(z) is noisy due to observational errors 
(mainly, due to variance in supernovae luminosity). There- 
fore, the smoothing of data over some interval Az is required 
(Az may depend on z). The value of Az is determined by 
estimated errors and by the required accuracy with which 
we want to determine H(z). Of course, the resulting H(z) 
will be smoothed, too, as compared to the genuine one. Note 
that our presentation here is very similar to that in Tegmark 
(2002). 

Instead of actually dividing a measured range of z 
into intervals, one may parametrize H(z) by some fitting 
curve which depends on a number of free parameters. This 
leads to model-independent parametric reconstruction of 
H(z), pde(z), w c h(z) and other quantities. It is clear that 
the number of free parameters N in such a fit just de- 
fines the equivalent smoothing interval Az (in particular, 
Az — z max /N if Az is chosen to be independent of z and 
we are considering the function H(z)/Ho, so that its value 
at z = is known exactly). Thus, the parametrization is 
equivalent to some kind of smoothing, with the actual way 
of smoothing (weighting) depending on the functional form 
of the parametric fit used. This refers even to such sophisti- 
cated methods as the 'principal-component' approach used 
in Huterer & Starkman (2002). Since decreasing Az (increas- 
ing N) results in a rapid growth of errors (AH(z) oc (5z)~ 3 '" 2 
directly follows from Eq. (34), c.f. Tegmark (2002)), for a 
given Zmax there is no sense in taking N to be large - this 
will merely result in the loss of accuracy of our reconstruc- 
tion. Thus, we will consider only 3-parametric fits for H(z) 
(these will correspond to 2-parametric fits for w(z)). 

After the discovery that the universe is accelerating, 
many different fitting function approaches were suggested 
and some are summarized below. 

• Polynomial Fit to Dark Energy : 

In this paper, we reconstruct dark energy using a very 
effective ansatz introduced in Sahni et al. (2003) in which 
the dark energy density is expressed as a truncated Taylor 
series polynomial in x — l + z, pr>E = A\ + A 2 x + A3X 2 . 
This leads to the following ansatz for the Hubble parameter 



which, when substituted in the expression for the luminosity 
distance (27), yields 



Dl 
l + z 



c 



±1 



VOml' 3 + A[ + A 2 X + A 3 X 2 



(36) 



The values of the parameters Ai,A 2 ,A3 are obtained 
by fitting (36) to supernova observations by means of a 
maximum likelihood analysis discussed in the next section. 
There are obvious advantages in choosing the ansatz (35) 
namely, it is exact for the cosmological constant w = — 1 
(A 2 — A3 = 0) as well as for quiessence with w = —2/3 
{Ai = A3 = 0) and w = -1/3 (Ai = A 2 = 0). Further- 
more, the presence of the term Q. m x 3 in (35) ensures that 
the ansatz correctly reproduces the matter dominated epoch 
at early times (z 3> 1). The presence of this term also allows 
us to incorporate information pertaining to the value of the 
matter density and, as we shall soon demonstrate, permits 
elaborate statistical analysis with the introduction of priors 
on Q m . 

The statefinder pair for the polynomial fit (35) can be 
written in terms of x — 1 + z as follows 



r(x) 
s(x) 



n m x 3 + A-! 



tt m x 3 + Ai + A 2 x + A3X 2 ' 

2(A 2 x + A 3 x 2 ) 
3(3^i + 2A 2 x + A3X 2 ) ' 



(37) 
(38) 



It is also straightforward to obtain expressions for the cos- 
mological parameters q and to by substituting (35) in (4) and 
(6) respectively. 

In figure 3 we show the maximum deviation between the 
exact value of the luminosity distance and the fit-estimated 
approximate value for a class of dark energy models. For 
LCDM (w — —1) and two quiessence models (w — —2/3, 
w — —1/3), the ansatz (36) returns exact values. (The ansatz 
is also exact for SCDM.) For the two tracker and Chaplygin 
gas models which we consider, the luminosity distance is 
determined to better than 1% accuracy for a conservative 
range in Q. m (0.2 < Q. m < 0.5). We therefore conclude that 
the polynomial fit (36) is very accurate and can safely be 
applied to reconstruct the properties of dark energy models. 

In this paper we will use the polynomial fit (36) to perform 
a model independent reconstruction of dark energy using 
the synthetic SNAP supernova data discussed earlier. Some 
details of our approach which involves the maximum likeli- 
hood method will be discussed in sections 4.2. Our results 
for the cosmological reconstruction of dark energy using the 
statefinder will be presented in section 5. 

Although we will mainly work with the polynomial ansatz 
(35) to reconstruct the properties of the statefinders, it is 
worthwhile to summarize some of the alternate approaches 
to the cosmological reconstruction problem. 

• Fitting functions to the luminosity distance Dl : 
An interesting complementary approach to the reconstruc- 
tion exercise is to find a suitable fitting function for the 
luminosity distance. Such an approach was advocated in 
Huterer & Turner (1999) and Saini et al. (2000). In Huterer 
& Turner (1999) a polynomial fit for the luminosity distance 
was suggested which had the form 



H{x) = H [Qmx 3 + A!+ A 2 x + A 3 x 2 ] " 



(35) 



Dl(z) 
l + z 



N 

E 

i=l 



aiZ 



(39) 
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Figure 3. The maximum deviation \ADj j /Dl\ between the 
actual value of the luminosity distance in the redshift range 
z = — 10 in a DE model and that calculated using the poly- 
nomial fit Eq (36). The solid line at AD^/D^ = represents 
models with w = — 1, w = —2/3, w = —1/3, for which the polyno- 
mial fit returns exact values . The dashed lines from top to bottom 
represent the tracker potential V{<j>) = Vo/4> a for a = 1 and 2 re- 
spectively. The dotted lines represent Chaplygin gas models with 
k = 0.5 and 2 (top to bottom). 



The ansatz (39) was examined in Weller & Albrecht (2002) 
who demonstrated that this approximation does not give 
an accurate reconstruction of the equation of state of dark 
energy. Similar conclusions will also be reached by us later 
in this paper in connection with the reconstruction of the 
statefinder pair using (39). 

A considerably more versatile and accurate fitting func- 
tion to the luminosity distance is (Saini et al. 2000) 



Dl 

x 



2 



■A ly /x- 1 + Ai 



A 2 x + Az^fx + 2 - Ai - A 2 - A 3 



x = 1+2 ,(40) 



where Ai, A2 and A3 are parameters whose values must be 
determined by fitting (40) to observations. Important prop- 
erties of this function are that it is valid for a wide range 
of models and that it exactly reproduces the results both for 
SCDM (fi m = 1) and the steady state model (Oa = 1). As 
demonstrated in Saini et al. (2000), an accurate analytical 
form for Dl allows us to reconstruct the Hubble parameter 
by means of the relation (34). Cosmological parameters in- 
cluding q(z), w(z), r(z), s(z) can now be easily reconstructed 
using (4), (6) and (32). 

• Fitting functions to the Equation of State : 
A somewhat different approach fits the equation of state 
of dark energy by the first few terms of a Taylor series ex- 
pansion (Weller & Albrecht 2002): 



wde(z) = 



WiZ 



(41) 



For N : 
D L 



1 the luminosity distance can be expressed as 



1 + 2 

Six 



c_ f 1+z 
Ho J, 



dx 



(1 



V^mX 3 + Q X 

Q m )x 3{1+wo) exp{3wi(a;- 1)} . 



(42) 



A modification of the above prescription was suggested in 
Gerke & Efstathiou (2002) which used a logarithmic expan- 
sion of the equation of state of dark energy: 



w(z) = too — a ln(l + 2) 



(43) 



where a — dw/d(lna). Yet another approach (Maor et al. 
2002) advocated a quadratic fit to the total equation of state: 

2 



Wt(z) = TOO + WlZ + W2Z 



(44) 



where the total equation of state, wt(z), is defined in terms 
of the equation of state of dark energy, to(z), as 



tot(z) 



w{z) 



1 + T%T ex P 



-3 



(45) 



Other approaches to the reconstruction problem can be 
found in Chiba & Nakamura (2000) ; Corasaniti & Copeland 
(2002); Linder (2003). 



4.2 Maximum Likelihood Estimation of 
cosmological parameters 

In order to determine how effective the statefinders are in 
discriminating between dark energy models, we adopt the 
method of maximum likelihood estimation to our recon- 
struction exercise. Supernova data is expected to improve 
greatly over the next few years. This improvement will be 
spurred by ongoing efforts by the Supernova Cosmology 
Project 1 and the High-z supernova search team 2 , as well 
by planned surveys such as the Nearby SN Factory 3 (300 
SNe at z < 0.1) and the SuperNova Acceleration Probe - 
SNAP 4 (~ 2000 SNe at z < 1.7). We shall use data sim- 
ulated according to the specifications of SNAP - a space 
based mission which is expected to greatly increase both 
the number of Type la SNe observed and the accuracy of 
SNe observations. 

SuperNova Acceleration Probe (SNAP) 
The SNAP mission is expected to observe about 2000 
Type la SNe each year, over a period of three years, ac- 
cording to the specifications given in Table 1. We assume 
a Gaussian distribution of uncertainties and an equidistant 
sampling of redshift in four redshift ranges. The errors in 
the redshift are of the order of Sz — 0.002. The statistical 
uncertainty in the magnitude of SNe is assumed to be con- 
stant over the redshift range ^ 2 ^ 1.7 and is given by 
""mag = 0.15. The systematic uncertainty limit is a 3ys = 0.02 
mag at redshift 2 = 1.5. For simplicity we assume a linear 
drift from a sys = at 2 = to a sys = 0.02 at 2 = 1.5, so 
that the systematic uncertainty on the model data is given 
by G 8ys {z) = (0.02/1.5)2. 

1 http://www-supcrnova.lbl.gov 

2 http://cfa-www.harvard.edu/cfa/oir/Research/supernova/HighZ.html 

3 http://snfactory.lbl.gov 

4 http://snfactory.lbl.gov 
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Optimizing the model with 2000 data points is some- 
what time consuming therefore we produced a smaller num- 
ber of binned SNe luminosity distances by binning the data 
in a redshift interval Az — 0.02. This interval is comparable 
to the statistical uncertainty in the redshift measurement 
of high-2 SNe due to the peculiar velocities of the galax- 
ies in which they reside, which is typically of the order of 
^peculiar ~ 1000 Kms -1 . In our experiment we smoothed 
the data in the first three redshift intervals in Table 1 by 
binning, the last interval had relatively fewer Sne and was 
left unbinned. The statistical error in magnitude, and hence 
in the luminosity distance is weighed down by the factor 
1/x/iVbm, where iVbin is the number of SNe in each bin. 

We use SNAP specifications to construct mock SNe cat- 
alogues. We may then use the method of maximum likeli- 
hood parameter estimation on this mock data to estimate 
the different cosmological parameters of interest. 

Maximum Likelihood Estimation: 

The observable quantity for a given supernova is its 
bolometric or 'apparent' magnitude m which is a measure of 
the light flux received by us from the supernova. To convert 
from m to cosmological distance, we use the well known 
relationship between the luminosity distance Dl and the 
bolometric magnitude 



m = M + 25 + 5\ogD L , 



(46) 



where Mo is the absolute magnitude of the SNe and the lu- 
minosity distance Dl is measured in the units of Mpc. (For 
Type la SNe, the typical apparent magnitude at z = 1 is 
about 25, which shows that we are dealing with very faint 
objects at that redshift.) Type la supernovae are excellent 
standard candles, and the dispersion in their apparent mag- 
nitude is CT mag = 0.15, which is nearly independent of the 
SN redshift. To relate this to the dispersion in the measured 
luminosity distance, we use Eq. (46) to obtain 



Cdist 

D L 



In 10 



= 0.069 



(47) 



While constructing mock SNe catalogues we shall assume 
that the errors in the luminosity distance are Gaussian with 
zero mean and dispersion given by the above expression (~ 
7%), the normalized likelihood function is therefore given by 



L( yi , Pk ) = n ( 



x exp 



27r<7dist(«») 
2 I &a\st{zi) 



(48) 



where the index i ranges from 1 to iVdat, which is the num- 
ber of supernovae in our sample, and we have denoted the 
fiducial supernova luminosity distance at a redshift z — Zi 
as yi = Dl(zi), where Dl(z) is the luminosity distance sim- 
ulated with SNAP specifications for a chosen background 
model using the (27). The p,'s are the parameters of the 
fitting function. (We shall mostly exploit the fitting func- 
tion (36) for which pi = Ai.) We maximize the Likelihood 
function L to obtain the Maximum Likelihood values of the 
parameters of the fitting function. In practice we minimize 
the negative of the log-likelihood, which is given by 



(49) 



where a constant term arising from the multiplicative factor 
is ignored. We are interested in estimating the statefinder 
pair r(z) and s(z) and the deceleration parameter q(z) from 
synthetic SNAP data. 

The priors that we have used for our reconstruction 
exercise are the following. 

The values of Ho and Mo (the absolute magnitude of 
SNe) are assumed to be known. We consider a flat uni- 
verse, so that the present day value of Qx is given by 
Qx = 1 — = Ai + Ai + A3. Also, when optimizing 
the model, we may assume priors on O m using information 
from other observations. This leaves only three free parame- 
ters (including Q m on which bounds can be specified). (Op- 
timizing without priors we found the variances of Ai to be 
much larger if no bounds were specified on f2 m .) 

Reconstruction of Cosmological Parameters 

Using the procedure described in detail above we now 
propose to reconstruct different cosmologically important 
quantities using SNAP data. We shall focus our atten- 
tion to the statefinder pair {r(z), s(z)}, the deceleration 
parameter q(z) and the cosmic equation of state w(z). 
Using SNAP specifications, we generated 1000 data sets 
{z^ ,D k Li , 0-^,0-% L .}, where the index k runs from 1-1000 
and the index i from 1 ~ 2000, with the LCDM as our 
fiducial model. For each of these experiments, the best-fit 
parameters, A 3 k (j — 1,2), for the polynomial fit to dark 
energy (35) were calculated. We then calculated r^{z) and 
s k (z) for each experiment from the calculated values A 3 k . 
The mean value of the statefinder pair and other cosmolog- 
ical quantities is computed as 



(r{z)) 



{*(*)) 



1UUU 

1=1 

1000 



(50) 



and so on for other quantities. Here the angular brackets 
denote ensemble average. We may also calculate the covari- 
ance matrix of these quantities at different redshifts which 
is given by 



Cr 



Cr 

C s 



(51) 



(52) 
(53) 
(54) 



where 

Crr = (r(z) 2 ) - (r(z)) 2 , 

C ss = (s(z) 2 ) - (s(z)f , 

C r . = (r(z)s(z))-(r(z))(s(z)) , 

and the angular averages are evaluates as in (50). 
5 RESULTS AND DISCUSSION 



From the results we can estimate the accuracy with which 
the ansatz recovers model independent values of different 
cosmological parameters, especially the statefinder pair in- 
troduced in Sahni et al. (2003). We can also determine 
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Table 1. Expectations from SNAP for a single year period of observation 



Rcdshift Interval 


z = 0-0.2 


z = 0.2-1.2 


z = 1.2-1.4 


z = 1.4-1.7 


Number of SNe 


50 


1800 


50 
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Figure 4. This figure shows 3<r confidence levels in the averaged statefinders (a) {s, f } and (b) {5, q}. The polynomial fit to dark energy, 
Eq (35) has been used to reconstruct the statefinders for an LCDM fiducial model with Q m = 0.3. The dashed line above the LCDM fixed 
point represents the family of quiessence models having w = constant. The dashed line below the LCDM fixed point shows Chaplygin 
gas models. It should be noted that the best-fit point in both panels (a) & (b) coincides with the LCDM fixed point (solid star). In the 
upper half of both panels, the solid rhombi correspond to tracker potentials V = Vo/4> a while triangles show ui=constant quiessence 
models. In the lower half of both panels, solid hexagons show Chaplygin gas models with different values of k. (The constant k gives the 
initial ratio between cold dark matter and the Chaplygin gas. Only models with k ^ f2 m /(l — Q m ) are permitted by theory, see Eq (15), 
(17).) All models, with the exception of the braneworld model, have Q m = 0.3 currently. The braneworld model is marked by a cross 
and corresponds to the DDG model (22) with Q m = 0.24 which best-fits current supernova data. Comparing the left and right panels 
we find that {s, q} is a slightly better diagnostic than {s, f } for tracker and quiessence models and can be used to rule out a constant 
equation of state ui ^ —0.9 at the 3<r level if the value of fi m is known exactly. 



whether this pair is useful in discriminating the cosmologi- 
cal constant model from other models of dark energy. 



5.1 Cosmological reconstruction for an LCDM 
fiducial model 

Synthetic supernova data are generated for a fiducial LCDM 
model with Qa = 0.7, Q m = 0.3 and assuming SNAP speci- 
fications summarized in the previous section. Next we deter- 
mine the statefinder pair and other cosmological parameters 
as functions of the redshift using the polynomial fit to dark 
energy (35). Our results can be represented in two comple- 
mentary ways. Firstly, we show the confidence levels in the 
r () — s space, where the subscript '0' denotes the present day 
value of the statefinders. We also find it useful to consider 



the integrated, averaged quantities: 



q(z) dz 



; max j q 
1 



r(z) dz 

1 

— / s(z)dz 

Zmax J 



(55) 
(56) 
(57) 



For the LCDM model, r and s do not evolve with time, 
therefore we find that f = 1 and s = 0. However for most 
other models of dark energy the statefinder pair evolves and 
the averaged quantities differ from their present day values. 
Due to averaging over redshift, the averaged parameters f, s 
are in many cases less noisy than ro, so. The maximum red- 
shift used for our reconstruction is z max = 1.7. One of the 
results of our analysis is that the deceleration parameter q 
is very well determined, see Fig. (9), therefore we also con- 
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Figure 5. This figure shows 3<r confidence levels in the statefindcrs: (a) {f, s}, (b) {q, s}, (c) {ro,so} and (d) {qo,so}- The fiducial 
model is assumed to be LCDM and, as in the previous figure, the polynomial fit to dark energy, Eq (35) is used to reconstruct the 
statcfindcr pairs. All notations are as in the previous figure. The current observational uncertainty in the value of the matter density is 
incorporated by marginalizing over the value of Q m - The dark grey outer contour shows results for the Gaussian prior f2 m = 0.3 ± Cfj m 
with OQ m = 0.05, the grey contour in the middle uses the Gaussian prior crQ m = 0.015, and the light grey contour is when Q m = 0.3 
exactly. Comparing panels (a) - (d) we find that {so, 90} is an excellent diagnostic of quintessence models which can be used to rule out 
a constant equation of state w —0.9 and tracker potentials V((p) oc <fi~ a , a J; 1, at the 3<x level even if Q m is known to an accuracy of 
only ~ 17%. It is important to note that of all statefinder pairs {so,<?o} is the least sensitive to the uncertainty in the value of f2 m . This 
is reflected in the fact that the 3<r confidence contour for {so,qo} with Q m = 0.3 ± 0.05 is not very much larger than the 3<x confidence 
level obtained if Q m is known exactly (Q m = 0.3). On the other hand the averaged statefinder pair {§, q} is a very good diagnostic of 
Chaplygin gas models and rules out models with n ^ 7 at the 3<x level if Q m = 0.3 ± 0.05. (The braneworld model marked by the cross 
can be ruled out by {so><?o} as well as {s, q}.) 
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Figure 6. Variation of < r(z) > with z for the cosmological 
constant model. Solid line shows best-fit < r(z) > averaged over 
all realizations calculated with the polynomial fit to dark energy, 
Eq (35), for the prior Q m = 0.3 exactly. The triple-dot-dashed 
line represents the exact value of < r >= 1 for the cosmological 
constant model. Shaded regions represent the Icr confidence levels 
for < r(z) >. The dark grey outer contour is for the Gaussian 
prior Q m = 0.3 ± cr fi m with o"n m = 0.05, the grey contour in 
the middle uses the Gaussian prior o-Q m = 0.015, and the light 
grey contour uses fl m = 0.3 exactly. The dotted, dashed and 
dot-dashed lines represent the exact values of r(z) for different 
constant w quicsscncc models, for kinessence models with the 
tracker potential V(<p) = Vo/4> a for different values of a, and for 
Chaplygin gas models with different n respectively. We see that 
all the model values plotted lie outside the la confidence level 
even for the most conservative prior of o"Q m = 0.05 at rcdshifts 
> 0.3. 



Figure 7. Variation of < s(z) > with z for the cosmological 
constant model. Solid line shows best-fit < s(z) > averaged over 
all realizations calculated with the polynomial fit to dark energy, 
Eq (35), for the prior Q m = 0.3 exactly. The triple-dot-dashcd 
line represents the exact value of < s >= for the cosmological 
constant model. Shaded regions represent the Icr confidence levels 
for < s(z) >. The dark grey outer contour is for the Gaussian 
prior Q m = 0.3 ± <xsi_ with crQ m = 0.05, the grey contour in 
the middle uses the Gaussian prior a"n TO = 0.015, and the light 
grey contour uses Cl m = 0.3 exactly. The dotted, dashed and 
dot-dashed lines represent the exact values of s(z) for different 
constant w quicsscncc models, for kinessence models with the 
tracker potential V(<j>) = Vo/4> a for different values of a, and for 
Chaplygin gas models with different re respectively. We see that 
all the model values plotted lie outside the Icr confidence level 
even for the most conservative prior of o"n m = 0.05. 



struct a second statefinder pair, {s,q}, which will be shown 
to be an excellent diagnostic of dark energy. 

Figure 4 shows the 99.73% confidence level in {s, f } (left 
panel) and {s,q} (right panel) for the fiducial LCDM model 
with O m = 0.3, Qa = 0.7. For comparison we also show 
values of f, s, q for quiessence, kinessence and Chaplygin gas 
models. From this figure we see that the statefinders can 
easily distinguish LCDM from: (i) quiessence with w >, —0.9 
(ii) the Chaplygin gas with k < 15 (iii) the quintessence 
potential V(<f>) oc 4>~ a ,a > 1 and (iv) the DDG braneworld 
models discussed in Deffayet et al. (2002). 

The above analysis assumed that the value of Q m is 
known exactly. However in practice it will be some time be- 
fore fl m is known to 100% accuracy and it is only natural 
to expect some amount of uncertainty in the observational 
value of this important physical parameter. We incorporate 
this uncertainty by marginalizing over the value of £l m . Two 
priors will be incorporated into our analysis, the weak Gaus- 
sian prior: Q. m = 0.3 ±0.05 and the stronger Gaussian prior: 
Q. m = 0.3 ±0.015. 



Figures 5 (a-d), show the confidence levels in the 
statefinder pairs {s, r}, {s, q}, {so,ro} and {so, go} respec- 
tively. For purposes of discrimination we also show the values 
of the r, s, q for quiessence, kinessence, Chaplygin gas and 
braneworld models. Figure 5(a) shows that the diagnostic 
{s, f } permits the LCDM model to be distinguished from 
quiessence models with w > —0.8, quintessence models with 
a > 1, Chaplygin gas models with k < 6 and braneworld 
models at the 99.73% confidence level and after applying 
the strong Gaussian prior of fl m = 0.3 ± 0.015. The dis- 
criminatory power of the statefinder clearly worsens for the 
weaker prior Q m = 0.3 ± 0.05. 

The situation can be dramatically improved if, instead 
of working with {s, r} we use the diagnostic {s, q} (see figure 
5 (b)). We find in this case that the fiducial LCDM model 
can be distinguished from quiessence with w > — 0.85 and 
the braneworld model at the 99.73% confidence level even 
for the weak prior Q. m = 0.3 ± 0.05. In Figures 5(c) & (d) 
we plot the confidence levels for current values of the pair 
{so,fo} and {so, go}. A few important points need to be 
noted here: 
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(a) (b) 

Figure 8. 3cr confidence levels in the parameter space f , s are shown for the cosmological constant model and the a = 1 tracker model 
using the polynomial fit to dark energy, Eq (35). The solid stars represent the model value of the parameter pair for the cosmological 
constant model and the a = 1 tracker kinessence models. The dashed line above the LCDM point represents the quiessence models, 
and that below the LCDM point represents different Chaplygin gas models. The solid triangles represent model values for constant w 
quiessence models, and the solid hexagons represent Chaplygin gas models with different values of k. Only those Chaplygin gas models 
with k ^ f! m /(l — Q m ) are allowed. For all the dark energy models, Q m = 0.3 is used. The ellipses represent the 3<r confidence levels 
in the f — s space for the exact prior Q m = 0.3. In (a), the dark grey contour in solid outline represents the confidence level for the 
cosmological constant fiducial model obtained when r, s are averaged over the redshift range z = to z = 1.7. The light grey contour 
with dotted outline is the confidence level for the a = 1 tracker kinessence model obtained when the averaging is over the entire redshift 
range. In (b), we show confidence levels for the cosmological constant model (dark grey contour) and for the a = 1 tracker model (light 
grey contour) with the averaging done from z = 1 to z = 1.7. Remarkably, using the statefinders {r , 5} one can rule out quintessence 
models with w > — 0.95 and Chaplygin gas models with k < 24 at 3<r if only very high redshift SNe belonging to the redshift interval 
1 Sj 2 ^ 1.7 are considered. The reason for this is that both r(z) and s(z) are determined to increasing accuracy at z > 1. Indeed, a 
'sweet spot' at z s ~ 1.4 ensures that both r(z s ) & s(z s ) are known with great accuracy at that point - see figures ( 6) & ( 7). 



(i) the semi-major axis of the confidence ellipse for {s, q} 
is tilted away from the dashed curve representing constant 
w models (quiessence). This enables the second statefinder 
pair {s, q} to be somewhat better at discriminating between 
LCDM and quintessence models than {r,s}. For instance, 
the current value {so,Qo} can discriminate the cosmological 
constant model from quiessence models having w >, — 0.9, 
kinessence models with a > 1, Chaplygin gas models with 
k < 2, and the braneworld model even after applying the 
weak Gaussian prior Q m — 0.3 ± 0.05. 

(ii) For Chaplygin gas models averaging over redshift 
considerably enhances the discriminatory prowess of both 
{s, f} and {s, q}. 

(iii) From figure (5d) we find that marginalization over 
Q m has only a small effect on the diagnostic {so, go} which 
contributes to making this statefinder pair a much better 
discriminator of dark energy than {r ,s } if the value of 
fi m is uncertain. 

Our results shown in figures 4 & 5 clearly demonstrate 
that both {r, s} as well as {s, q} are excellent diagnostics of 
dark energy with the latter being somewhat more sensitive 
than the former. 



We now proceed to examine the information content in 
the cosmological parameters when examined individually. 
In Figure 6, we plot the variation of the ensemble averaged 
value (r(z)) with redshift. The la error bounds are shown 
for two different priors on O m and for the case when the 
value of Q, m is known exactly. This figure shows that r(z) is 
a good diagnostic of dark energy and allows us to discrimi- 
nate (at the 68.3% CL) between different dark energy mod- 
els and LCDM. Discrimination improves at higher redshifts 
(z > 0.8) especially if the uncertainty in the value of O m is 
small. The sweet spot for this parameter, i.e. the point at 
which r(z) is most accurately determined, is at z ~ 1.4. (For 
earlier work on the sweet spot see Weller & Albrecht 2002; 
Huterer & Turner 1999; Huterer & Starkman 2002). 

Figure 7 shows the variation of the ensemble averaged 
value of the second statefinder (s(z)) with redshift. Again la 
errors for the two priors fl m = 0.3 ± 0.05, Q m = 0.3 ± 0.015 
and when the value of Q. m is known exactly (fi m = 0.3) are 
shown. We see that s(z) is determined even more accurately 
than r(z), and therefore can serve as a better diagnostic of 
dark energy. For the strong Gaussian prior Q. m = 0.3±0.015, 
(or when fi m is known exactly) the value of s is very well 
determined even at higher redshifts, its sweet spot being 
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Figure 9. Variation of < q(z) > with z for the cosmological 
constant model. Solid line shows best-fit < q(z) > averaged over 
all realizations calculated with the polynomial fit to dark energy, 
Eq (35), for the prior Q m = 0.3 exactly. Shaded regions repre- 
sent the la confidence levels for < q >. We find here that the 
use of exact Q m = 0.3 and the Gaussian prior f2 m = 0.3 ± crn m 
with o"Q m = 0.015 gives us almost the same bounds, represented 
by the light grey contour, the dark grey outer contour uses the 
Gaussian prior aQ m = 0.05. The dotted, dashed and dot-dashed 
lines represent the model values of < q > for different constant w 
quicssence models, for kinessence models with the tracker poten- 
tial V(4>) = Vo/<f> a for different values of a, and for Chaplygin gas 
models with different n respectively. The horizontal thick line rep- 
resents the accuracy with which the acceleration epoch is deter- 
mined in this fit. This figure demonstrates that q can be used as a 
discriminator between dark energy models at low rcdshifts < 0.5. 
Indeed the location of a sweet spot at z B ~ 0.2 demonstrates that 
the value of q{z s ) is known to remarkably good accuracy ! 

at z ~ 1.4. For the weak prior O m = 0.3 ± 0.05, s is not 
so well determined at high redshifts, but it is still accu- 
rate enough to distinguish between rival models of dark en- 
ergy. Two points are of interest here. Firstly, r and s are 
both much more accurately determined at higher redshifts 
if the value of Q. m is accurately known. This explains why 
the parameters {f, s} perform better as discriminators than 
{ro, so}- Secondly, the sweet spot for both these parameters 
appears at z ~ 1.4, only if the value of Q m is accurately 
known. Upon marginalizing over O m the sweet spot disap- 
pears both in the case of r(z) as well as in the case of s(z). 
Another point worth mentioning is that Chaplygin gas mod- 
els are much easier to rule out at high z than at low z, using 
either r(z) or s(z). As an illustration, neither ro nor so can 
distinguish a k — 2 Chaplygin gas model from LCDM (with 
identical tl m ) at the la level. However the averaged-over- 
redshift statefinders f , s can do so quite easily even at the 
3cr level, as demonstrated in figures 4 and 5. 

Figure 8 shows how sweet spot information can be used 
to improve the statefinders as a diagnostic tool. For both r 



and s the sweet spot appears at high redshifts. Therefore, 
one expects that the discriminatory prowess of the statefind- 
ers will improve considerably if only data at z ^ 1 is con- 
sidered. This is indeed the case. Figure 8 shows 3a confi- 
dence levels in {f, s} for two cases: (a) the statefinder pair 
is averaged over the full redshift range < z < 1.7, (b) 
the statefinder pair is averaged over the high redshift range 
1 < z < 1.7; Q m = 0.3 for both cases. The dark grey el- 
lipses represent the confidence level for the LCDM model, 
and the light grey ellipses represent the confidence level for 
the a = 1 kinessence model. We see that there is a dramatic 
improvement in the determination of the statefinder pair in 
figure 8(b) where the statefinders have been averaged only 
for 2 ^ 1. From figure 8 (a) we see that {f , s} can discrimi- 
nate between LCDM and quiessence models with w >, —0.90 
and Chaplygin gas models with k < 15, whereas 8(b) shows 
that {f, s} can discriminate between LCDM and quiessence 
models with w _> — 0.95 and Chaplygin gas models with 
k < 24 ! We therefore conclude that high redshift super- 
novae can play a crucial role in constraining properties of 
dark energy and our results support the views expressed in 
Linder & Huterer (2003). We must however note that in 
order to use sweet spot information optimally the value of 
fi m must be known to very high accuracy. Indeed, for the 
Gaussian prior of O m = 0.3 ± 0.05, a consideration of only 
high redshift supernovae does not lead to any improvement 
in the results. This is because, as seen from Figures 6 & 
7, after marginalization over O m , the sweet spot for both 
r(z) and s(z) disappears. The second point to note is that 
the angle of inclination of the semi-major axis of the ellipse 
to the w = constant curve (quiessence) appears to depend 
upon the redshift range over which the statefinder pair is 
being averaged. 

Figure 9 shows the variation of the mean deceleration 
parameter q(z) with redshift. We see that q(z) is very well 
determined over the entire range ^ z ^ 1.7. This justifies 
our choice of {s, q} as the second statefinder pair. Indeed, 
the behaviour of r(z) and s(z) on the one hand and q(z) on 
the other, is in some ways complementary. While both r(z) 
and s(z) are determined to increasing accuracy at higher red- 
shifts, the deceleration parameter is very well determined at 
lower redshifts and the sweet spot for this parameter appears 
at the redshift z s ~ 0.25. It is interesting that, in sharp con- 
trast to what was earlier observed for r and s, the sweet spot 
in q(z) persists even after we marginalize over Q, m ! From 
figure 9 we can also determine the value of the acceleration 
epoch (the redshift at which the universe began accelerat- 
ing). We find that the acceleration epoch is determined quite 
accurately: z(q = 0) = 0.66 ± 0.06. 

Figure 10 shows maximum likelihood contours for the 
pair {s, q} where q has been averaged over the redshift inter- 
val < z < 0.4 while s has been averaged over 1 < z < 1.7. 
This figure clearly demonstrates that {s, q} is an excellent 
diagnostic of dark energy since it can distinguish LCDM 
from quiessence models with w ^ —0.95 on the one hand 
and from Chaplygin gas models with k ^ 25 on the other 
(both at the 99.73% CL). Figures 8 and 10 show that the 
ability of the averaged statefinder pairs {s, f } and {s, q} to 
discriminate between dark energy models is comparable if 
the value of Q. m is known exactly. (As demonstrated earlier 
in figure (5), {s, q} is a more sensitive diagnostic than {r, s} 
if we marginalize over Q m -) 
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Figure 10. 3tr confidence levels in the parameter space q, s are shown for the cosmological constant model and the a = 1 tracker model 
using the polynomial fit to dark energy, Eq (35). The solid stars represent the model value of the parameter pair for the cosmological 
constant model and the a = 1 tracker kinessence models. The dashed line above the LCDM point represents the quiessence models, 
and that below the LCDM point represents different Chaplygin gas models. The solid triangles represent model values for constant w 
quiessence models, and the solid hexagons represent Chaplygin gas models with different values of k. Only those Chaplygin gas models 
with k ^ f! m /(l — Q m ) are allowed. For all the dark energy models, Q m = 0.3 is used. The ellipses represent the 3<r confidence levels 
in the q — s space for the exact prior Q m = 0.3. fn (a), the dark grey contour in solid outline represents the confidence level for the 
cosmological constant fiducial model obtained when q,s are averaged in the redshift range z = to z = 1.7. The light grey contour with 
dotted outline is the confidence level for the a = 1 tracker kinessence model obtained when the averaging is over the entire redshift 
range, fn (b), we show confidence levels for the cosmological constant model (dark grey contour) and for the a = 1 tracker model (light 
grey contour) with the averaging done from z = to z = 0.4 for q and from z = 1 to z = 1.7 for s. The reason for choosing these ranges 
is that q is extremely well-determined at low redshifts, with a sweet spot at z s ~ 0.25, and s is accurately determined at high redshifts 
with a sweet spot at 2 S ~ 1.4 — see figures ( 9) & ( 7). Using the second statcfinder pair {q, s} one can rule out quintessence models with 
w > — 0.95 and Chaplygin gas models with k < 25 at 3<r if only very high redshift SNe belonging to the redshift interval 1 $J z ^1.7 
are considered for s and low redshift SNe in the interval ^ z ^ 0.4 arc considered for q. 



Figure 11 shows how the equation of state w(z) varies 
with redshift. We see that, although the equation of state is 
determined remarkably well at small redshifts, cosmological 
reconstruction of w(z) steadily worsens with z and deteri- 
orates rapidly beyond z ~ 1. With the most conservative 
prior f2 m = 0.3 ± 0.05, the ansatz (35) can distinguish be- 
tween the cosmological constant model and the quiessence 
model with w = —0.96 at the 3a level provided we restrict 
ourselves to low redshifts z < 0.4. For higher redshifts, the 
LCDM model cannot be distinguished from the w = —0.8 
model beyond z ~ 1.1 (after marginalizing over Q m with the 
prior fi m = 0.3 ±0.05). In the ideal case when O m = 0.3 ex- 
actly, LCDM and the w — —0.8 model can be distinguished 
up to z ~ 1.3 but not beyond. (In this case the best-fit w(z) 
is biased beyond z ~ 0.7 and takes on a lower value than 
the fiducial w — —1.) Somewhat surprisingly, although w(z) 
and q(z) are related through (6) and therefore carry essen- 
tially the same information, even a cursory examination of 
figures 9 and 11 reveals that the ansatz (35) does not re- 
construct w(z) to the same accuracy as it reconstructs q(z). 
However, like q(z), w(z) is reasonably well determined at low 
redshifts, having a sweet spot at z ~ 0.25. The sweet spot 



persists when we marginalize over Q m using the Gaussian 
prior Q. m = 0.3 ±0.015 but disappears when the uncertainty 
in O m is increased using the prior Q m = 0.3 ± 0.05. 



5.2 Cosmological reconstruction for a tracker 
model 

We now briefly examine the accuracy of the statefinder pair 
and the ansatz (35) in determining the statefinder pair for a 
fiducial model other than LCDM. We know that the ansatz 
returns exact values for the cosmological constant w = — 1 
as well as for quiessence having the constant equation of 
state w = —2/3 and w = —1/3 (see figure 3). It is therefore 
important to study the accuracy of the statefinder in recon- 
structing the properties of dark energy in models in which 
both the dark energy density as well as the equation of state 
vary with time and for which the ansatz (35) is approximate. 
For this purpose we shall work with a fiducial dark energy 
model which evolves under the influence of the tracker po- 
tential V — Vo/4> and use the ansatz (35) in tandem with 



Exploring the Expanding Universe and Dark Energy using the Statefinder Diagnostic 17 



A 

s 

V 




w = 0.0 




QUINTESSENCE 



> ot~ 1 



LCDM 



/c=6« 



CHAPLYGIN GAS 



K = 2» 

I 

/ 

K=1 » 



*=3,{7,,-- 



0.5 



Figure 11. Variation of < mj(z) > with z for the cosmological 
constant model. Solid line shows best-fit < w(z) > averaged over 
all realizations calculated with the polynomial fit to dark energy, 
Eq (35), for the prior Q m = 0.3 exactly. The dot-dashed line 
represents the exact value of < w >= —1 for the cosmological 
constant model. Shaded regions represent the Icr confidence lev- 
els for < w >. The dark grey outer contour is for the Gaussian 
prior Q m = 0.3 ± o~q with o~n = 0.05, the grey contour in 
the middle uses the Gaussian prior cr(i m = 0.015, and the light 
grey contour uses Q m = 0.3 exactly. The dotted and dashed lines 
represent the model values of < w > for different constant w 
quicssence models and for kinessence models with the tracker po- 
tential V{4>) = Vo/(j> OL for different values of a respectively. We see 
that w can distinguish between the cosmological constant model 
and other dark energy models only at z < 1. 



the statefinders (32) to reconstruct the properties of dark 
energy. 

Figure 12 shows our results in terms of 3a confidence 
levels in {f, s}. We find that, if the value of f2 m is known 
to reasonable accuracy (Q, m — 0.3 ± 0.015) then the aver- 
aged statefinder pair {f, s} is able to distinguish the model 
V = Vo/4> from the model V = Vo/(f> 2 at the 3ct level. 
As expected, a large uncertainty in the current value of 
the matter density reduces the efficiency of this diagnos- 
tic and the two models V = Vo/cp and V = Vo/4> 2 cannot be 
clearly distinguished if the uncertainty in Q m is increased to 
n m = 0.3 ±0.05. 

Figures 13 & 14 show the performance of the ensemble 
averaged statefinders (r(z)) and (s(z)). As in the earlier case 
when our fiducial model was assumed to be LCDM, we find 
that the errors in r(z) and s(z) are small. However a slight 
bias in the determination of the statefinders exists at low 
redshifts so that for z < 0.4 the value of the best fit (r) ((s)) 
is larger (smaller) than the exact fiducial value. Averaging 
over the entire redshift range can significantly reduce this 
bias and we conclude that the ansatz (35) works well even 
for those dark energy models for which it docs not return 



r 

Figure 12. 3cr confidence levels in the parameter space f, s are 
shown for the kinessence model with tracker potential V = Vo /4> a 
for a = 1, for the polynomial fit to dark energy, Eq (35). The 
solid star represents the model value of the parameter pair for the 
cosmological constant model. The dashed line above the LCDM 
point represents the quiessence models, and that below the LCDM 
point represents different Chaplygin gas models. The solid rhombi 
represent tracker kinessence models with different a, the solid tri- 
angles represent the constant w quicssence models, and the solid 
hexagons represent Chaplygin gas models with different values 
of k. Only those Chaplygin gas models with k ^ S7 m /(1 — Q m ) 
are allowed. For all the dark energy models, Q m = 0.3 is used. 
The square represents the best-fit for the a = 1 fiducial model. 
The shaded ellipses represent the 3cr confidence levels in the f — s 
space. The dark grey outer contour is for the Gaussian prior with 
dfi m = 0.05, the grey contour in the middle uses the Gaussian 
prior <rn m = 0.015, and the light grey contour uses Q m = 0.3 
exactly. 



exact values. One should also note the reappearance of the 
sweet spot for the statefinders r & s in the figures 13 & 
14. For the statefinder r(z) the sweet spot appears at z ~ 
1.6, while for s(z) the sweet spot is at z ~ 1.2. As in the 
case of the LCDM model, one can try and constrain the 
properties of dark energy further by evaluating the averaged 
statefinders using SNe data only from z > 1. Our results 
shown in figure 8 demonstrate that the confidence ellipse 
for {s, f } becomes much smaller when the averaging is done 
over the redshift range 1 ^ z < 1.7 than when the averaging 
is over the entire redshift range. 

The performance of the deceleration parameter (q{z)) 
for this quintessence model is shown in the figure 15. We 
see that the deceleration parameter is very accurately de- 
termined. The sweet spot for the deceleration parameter oc- 
curs at lower redshifts, at z ~ 0.2, and by combining higher 
redshift data in determining s with lower redshift data for 
determining q we can significantly improve the errors on the 
second statefinder pair {s,q}, as demonstrated in figure 10. 
As was the case for the LCDM model, the sweet spot gradu- 
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Figure 13. Variation of < r(z) > with z for the kincsscncc model 
with the tracker potential V(</>) = Vo/<f> a for a = 1. Solid line 
shows best-fit < r(z) > averaged over all realizations calculated 
with the polynomial fit to dark energy, Eq (35), for the prior 
f^m = 0.3 exactly. Shaded regions represent the Icr confidence 
levels for < r >. The dark grey outer contour is for the Gaussian 
prior Q m = 0.3 ± ctq with o"n m = 0.05, the grey contour in 
the middle uses the Gaussian prior cn m = 0.015, and the light 
grey contour uses Q m = 0.3 exactly. The dotted, dashed and dot- 
dashed lines represent the model values of < r > for different 
constant w quiessence models, for tracker kinessence models for 
different values of a, and for Chaplygin gas models with different 
k respectively. The thick solid line represents LCDM. 



Figure 14. Variation of < s(z) > with z for the kinessence model 
with the tracker potential V(</>) = Vo/4> a for a = 1. Solid line 
shows best-fit < s(z) > averaged over all realizations calculated 
with the polynomial fit to dark energy, Eq (35), for the prior 
f^m = 0.3 exactly. Shaded regions represent the la confidence 
levels for < s >. The dark grey outer contour is for the Gaussian 
prior Q m = 0.3 ± o"r> m with a-Q m = 0.05, the grey contour in 
the middle uses the Gaussian prior crj m = 0.015, and the light 
grey contour uses f2 m = 0.3 exactly. The dotted, dashed and dot- 
dashed lines represent the model values of < s > for different 
constant w quiessence models, for tracker kinessence models for 
different values of a, and for Chaplygin gas models with different 
k respectively. The thick solid line represents LCDM. 



ally disappears if we incorporate the prevailing uncertainty 
in the value of the matter density by marginalizing over large 
values of fi m . 



5.3 Cosmological reconstruction using other 
fitting functions 

For comparison, we also carry out the reconstruction exer- 
cise using two of the fitting functions described in section 
4.1. In figure 16 we show the results for {r,s} using the 
polynomial fit to the luminosity distance (39) with N = 5. 
In this case, because of the nature of the ansatz, it is not 
possible to place any priors on fi m . We find that this ansatz 
does not perform well for the statefinder pair. Firstly, it does 
not determine f , s with the accuracy seen in the case of the 
polynomial fit to dark energy. Secondly, the best-fit value 
for {f, s} is biased with respect to the fiducial LCDM value. 
Additionally, the errors on both r and s are unacceptably 
large due to which one cannot distinguish between the cos- 
mological constant model and kinessence models with a < 6, 
quiessence models with w < — 0.4, and Chaplygin gas mod- 
els with k > 2at the 3cr confidence level. Even at la (68.3% 
CL), one can only discriminate LCDM from kinessence mod- 
els with q ^ 3, quiessence models with w > — 0.6, and 



Chaplygin gas models with n < 3. We therefore conclude 
that the polynomial fit to the luminosity distance (39) is 
not very useful for the reconstruction of the statefinders. 

We also carry out a similar reconstruction exercise using 
the polynomial fit to the equation of state (41) with N = 1. 
This ansatz can accommodate priors on Q m and we expect it 
to perform better than the polynomial fit to the luminosity 
distance. Indeed, figure 17 clearly demonstrates that a two 
parameter Taylor expansion in the equation of state is better 
than a five parameter expansion in the luminosity distance 
(Our results in this context support the earlier observations 
of Weller & Albrecht 2002). From figure 17 we find that the 
ansatz (41) can distinguish the cosmological constant model 
from quiessence models with w > — 0.6, kinessence models 
with a ^ 3, and Chaplygin gas models with k < 3 at the 
99.73% confidence limit after we have marginalized over Q. m 
with a Gaussian prior of fl m = 0.3 ± 0.05. However a com- 
parison of figure 17 with figure 5 shows that the equation of 
state expansion (41) is not quite as accurate as the polyno- 
mial fit to dark energy (35) in reconstructing the statefinder 
pair {f, s}. We therefore conclude that the statefinders can 
be reconstructed using several complementary methods. The 
polynomial fit for dark energy (35), by providing a good re- 
construction of the parameters f , s, q, can successfully be 
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Figure 15. Variation of < q(z) > with z for the kinessence model 
with the tracker potential V(<j>) = Vo/<f> a for a = 1. Solid line 
shows best-fit < q(z) > averaged over all realizations calculated 
with the polynomial fit to dark energy, Eq (35), for the prior 
fim = 0.3 exactly. Shaded regions represent the la confidence 
levels for < q >. We find here that the use of exact f2 m = 0.3 and 
the Gaussian prior Q m = 0.3 ± o"r> m with aQ m = 0.015 gives us 
almost the same bounds, represented by the light grey contour, 
the dark grey outer contour uses the Gaussian prior aQ m = 0.05. 
The dashed and dot-dashed lines represent the model values of 
< s > for tracker kinessence models for different values of a, and 
for Chaplygin gas models with different k respectively. The thick 
solid line represents LCDM. 

used for the model independent reconstruction of dark en- 
ergy. 



6 CONCLUSIONS AND DISCUSSION 

This paper contains an in depth study of properties of the 
statefinder diagnostic introduced in Sahni et al. (2003). The 
statefinder pairs {r, s} and {s, q} have the potential to suc- 
cessfully discriminate between a wide variety of dark energy 
models including the cosmological constant, quintessence, 
the Chaplygin gas and braneworld models. The statefind- 
ers play a particularly important role in the case of modi- 
fied gravity theories such as string/M-theory inspired scalar- 
tensor models and braneworld models of dark energy, for 
which the equation of state is not a fundamental physical 
entity and therefore does not provide us with an adequate 
description of the accelerating universe. Our results, sum- 
marized in figures 1 and 2, show that the statefinders r, 
s considerably extend and supplement traditional measures 
of cosmological dynamics such as the deceleration parameter 
q. To give a concrete example of how this can happen con- 
sider two (or more) cosmological dark energy models which 
have identical (hence degenerate) values of q . Although 
such models will have the same current value of a/ a, the 
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Figure 16. 3a confidence levels in the parameter space f, s are 
shown for the cosmological constant model using the polyno- 
mial fit to luminosity distance, Eq (39). The solid star represents 
the model value of the parameter pair for the cosmological con- 
stant model. The dashed line above the LCDM point represents 
the quiessence models, and that below the LCDM point repre- 
sents different Chaplygin gas models. The solid rhombi represent 
tracker kinessence models with different a, the solid triangles rep- 
resent model values for constant w quiessence models, and the 
solid hexagons represent Chaplygin gas models with different val- 
ues of k. Only those Chaplygin gas models with k ^ f2 m /(l — f2 m ) 
are allowed. For all the dark energy models, Q m = 0.3 is used. 
The best-fit point for the reconstruction is represented by the solid 
square. The ellipses represent the la, 2a, 3a confidence levels in 
the f — s space. 

value of the third derivative a (hence r & s) will in gen- 
eral be different in both models. The statefinder pairs {r, q} 
and {s, q} therefore provide us with a 'phase-space' picture 
of dark energy which distinguishes dynamical dark energy 
models both from each other as well as from the cosmological 
constant and helps break cosmological degeneracies present 
in rival models of dark energy. The statefinder s is remark- 
ably sensitive to the total pressure of all forms of matter 
and radiation in the universe. As a result s remains sensi- 
tive to the presence of dark energy even at moderately high 
redshifts z ~ 10, when the universe is matter dominated. 

Forthcoming space-based missions such as SNAP are 
expected to greatly increase and improve the current Type 
la supernova inventory. Anticipating this development we 
have carried out a maximum likelihood analysis which com- 
bines the statefinder diagnostic with realistic expectations 
from the SNAP experiment. Our results, summarized in 
figures 4 - 10, show that both {r, s} as well as {s, q} are 
excellent diagnostics of dark energy. If the value of Q m is 
known exactly, then the averaged-over-redshift statefinder 
pair {s, q, } can distinguish between the cosmological con- 
stant model (w = —1) and a dark energy model having 
w = —0.9 at the 99.73% CL. It can also distinguish (at the 
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are better constrained at lower z (z < 0.4). One is there- 
fore tempted to construct a new diagnostic {s, q}, where s 
is averaged over the redshift range 1 ^ z ^ 1.7 whereas <? 
is averaged over the redshift range < z < 0.4. From fig- 
ures 8 & 10 we find that {s, q} provides an extremely potent 
diagnostic of dark energy since it can distinguish a fiducial 
LCDM model from a dark energy model with w > —0.95 on 
the one hand, and from a Chaplygin gas model with k < 25 
on the other, at the 99.73% confidence level. 

We therefore believe we have convincingly demon- 
strated that the statefinder pair {s, f } and {s, q} provide 
an excellent diagnostic of dark energy which will be used to 
successfully differentiate between the cosmological constant 
and dynamical models of dark energy. 
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Figure 17. 3c confidence levels in the parameter space f,s are 
shown for the cosmological constant model using the polyno- 
mial fit to equation of state, Eq (41). The solid star represents 
the model value of the parameter pair for the cosmological con- 
stant model. The dashed line above the LCDM point represents 
the quiessence models, and that below the LCDM point repre- 
sents different Chaplygin gas models. The solid rhombi represent 
tracker kinessence models with different a, the solid triangles rep- 
resent model values for constant w quiessence models, and the 
solid hexagons represent Chaplygin gas models with different val- 
ues of re. Only those Chaplygin gas models with re ^ Q m /(l — Cl m ) 
are allowed. For all the dark energy models, Q m = 0.3 is used. 
The best-fit point lies on the LCDM point for this reconstruc- 
tion. The shaded ellipses represent the 3cr confidence levels in the 
f — s space. The dark grey outer contour is for the Gaussian prior 
n m = 0.3±o-Q m with <7o. m = 0.05, the grey contour in the middle 
uses the Gaussian prior o"n m = 0.015, and the light grey contour 
uses il m = 0.3 exactly. 



same level of confidence) the cosmological constant (LCDM) 
from Chaplygin gas models with k < 7. 

Keeping in mind the fact that the current observational 
data determine Q. m to a finite level of accuracy, we have 
probed how well the statefinder fares as a diagnostic after 
one incorporates the prevailing uncertainty in the value of 
the matter density by marginalizing over values of £l m which 
are uncertain. Somewhat surprisingly, the statefinder fares 
rather well even for the relatively weak prior Q m = 0.3±0.05. 
In this case, by employing the diagnostic {so, go}, the LCDM 
model can be differentiated from the w = —0.9 model on the 
one hand, and from the tracker potential V((f>) oc (p^ 1 and 
the DDG braneworld model (Deffayet et al. 2002) on the 
other, at the 99.73% CL. 

Finally we should mention that of the two statefinders 
s(z) appears to be better constrained by observations espe- 
cially if the uncertainty in 57 m is small. Interestingly both 
r(z) and s(z) show less scatter at higher redshifts (z > 1) 
and thereby complement the behaviour of the deceleration 
parameter q(z) and the cosmic equation of state w(z) which 
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